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Abstract 

If calculated in the standard way, the cross section for the collision of two 
unstable particles turns out to diverge. This is because part of such a 
cross section is proportional to the size of the colliding beams. The effect 
is called the "linear beam size effect" . We present a way of including this 
linear beam size effect in the usual Monte Carlo integration procedure. 
Furthermore we discuss the gauge breaking that this may cause. 



1 Introduction 



The cross section for the collision of two unstable particle generally diverges. 
This happens for instance in the Feynman graph 




W+(q 3 ) 



(1) 



The lower half of this graph looks like the decay of a muon. Consequently the 
kinematics of the process allows the momentum k to be on its mass shell. After 
all that is what one gets from the decay of a muon: a muon neutrino on its 
mass shell. The factor l/(fc 2 + ie) that occurs in the matrix element causes a 
divergence of the total cross section. 

In pP and 0] this problem has been studied in detail, and it has been shown 
that this divergence is softened into a finite peak if the incoming particles are 
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described carefully enough. In this paper we give a prescription for including 
this peak in Monte Carlo simulations. Typically such modifications may result 
in a violation of gauge invariance in the amplitude. We study this effect in 
detail. 



2 Asymptotic States 

In the context of scalar particles, Veltman [5] has shown that the S-matrix 
satisfies unitarity and causality only if one restricts the in/out states to the sta- 
ble particles. Because of this, when considering collisions of unstable particles, 
we should use Feynman graphs that take the production process of the unsta- 
ble particle into account. We are going to show that we actually do not need 
to worry about this as long as the wave packets of the unstable particles are 
much smaller in size than the decay length of the unstable particle. A complete 
amplitude for the production and collision of two muons looks like 

A=(2n)H f ffi ^ e -m^/% e -^W%f«(p> +p > 2 _ gl _ g2 ) 



i> pt (pi) = / dpu dp' b <b Pa ( P ' a )cb Pb {p' b ){2itf5\p' a +p' b - p[)M pmd (3) 



where 



may be viewed as the wave function of the unstable particle. Notations like 
4>p a (p' a ) stand for the wave function of a particle that is peaked in momentum 
space around the value p a evaluated at p' a The above expression for the wave 
function of an unstable particles assumes that the unstable particle is produced 
in a two-to-one process. We assume this only for the sake of simplicity of 
notations. If there are other outgoing or incoming particles their wave functions 
can easily be added. Also note the factors e~ lTiPi ' Pi ^ m ^ . These factors are 
translations of the wave function. The point of these translations is that they 
make sure that the unstable particles are produced away from the spot where 
they collide. The invariant distance that the unstable particle travels before 
colliding is tj. 

Now we are going to use the assumption that the wave packets are much 
smaller than the decay length. This has as a consequence that in momentum 
space the wave packets are much broader than the decay length. Because of 
this we may assume that they are constant functions of (pi) 2 resp (p'2) 2 over a 
range of several times m^T^. Therefore it is possible to integrate the expression 
for A given above over the values of (p[) 2 and {p'2) 2 ■ We only have to integrate 
the factors contained in the quantity F that is defined to be given by 

p = e -«i"ipi-pi/«V e -iT2P2-p 2 /nV 



[(p' 2 ) 2 - ml + im^Tf,} [(p[) 2 - + im^T^] ' 

(4) 
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We integrate along a path parameterized as p[ 2 {t). This parameterization is 
done according to 

p[(t) =pi(0)+ic; 

k'(t) =k'(0)+tc; (5) 
P2(t) =p' 2 (0)-tc. 

If we want to integrate over the value of (pi) 2 we choose c to be a four vector 
that is a linear combination of p\, p 2 and k such that it is orthogonal to the 
latter two vectors but not to pi . This parameterization is chosen such in order 
to satisfy momentum conservation and furthermore to be on a constant fc 2 -plane 
in order not to get difficulties with the singularity at k 2 = 0. After doing this 
integral and an analogous one over the value of (p' 2 ) 2 , we find the result 

A=(2n) i i { dp[dp 2 e~ iTlp i- pi/m »e~ r » Tl/2 e- iT2p 2- p2/m »e- r » T2/2 

( 6 ) 

<5 4 (p'i +P 2 - qi - q2)^p 2 (p 2 )u l (p' 2 )u i (p 2 )M co \iu j (p[)u 3 (p[)^ pi (p' 1 ). 

This is (except for the decay factors e~ TflTi ^ 2 ) exactly the same as if we had 
started with incoming muons on their mass shell. The conclusion is that if we 
have wave packets that are much smaller than the decay length of the unstable 
particles we may treat them as if they were asymptotic states. 

This has no bearing on the question of gauge invariance. Matrix elements 
depend on the masses of the particles. If masses are chosen such that the muon 
is no longer unstable (by assuming the electron to be heavier than the muon, so 
that the decay is forbidden), the matrix element is gauge invariant, so it must 
also be if masses are chosen accordingly to their measured values. 



3 The Linear Beam Size Effect 

We observed that the divergence in the cross section is caused by a peak in 
the matrix element in momentum space. A sharp peak in momentum space 
means a long range effect in position space. Indeed, the decay product of a 
decaying muon can travel over arbitrary distances. The consequence is that 
the cross section becomes proportional to the size of the beam. In colliders the 
longitudinal beam size is much larger than the transverse one. Consequently, the 
cross section is actually proportional to the transverse beam size, to be denoted 
by a. The precise definition of this quantity can be found in [I]. In the same 
paper a more rigorous version of this qualitative argument was given. In 0] it 
was shown that the quantities used in the rigorous argument can replaced by 
covariant ones. 

The part of the cross section proportional to the beam size is given by 

CTscmi-sing = CUT f da IC( i— :S(k 2 — TO 2 ), (7) 

J \k±\ 
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where tr re d is the cross section with the offending propagators [k 2 — m? ± ie) _1 
removed. k± is by definition given by k + api + 0p2 with a and (3 chosen such 
that k± ■ pi.2 = 0. 

The above formula gives the part of the cross section proportional to the 
beam size, but it would be more convenient if the linear beam size effect could 
be incorporated in the usual Monte Carlo integration procedure. This can indeed 
be done by doing the substitution 



k 2 — m 2 ±ie k 2 — m 2 ± i\k±\/a 
If we use the approximation 

(fc 2 - to 2 ) 2 + \k ± \ 2 /a 2 ~ \k~Z\ ^ ~ m2) ' (9) 

these two expressions become equal. This approximation only needs to be valid 
around the peak at k 2 = m 2 , which will generally be the case. The only property 
that is needed for this to be true is that the reduced cross section da Te d does 
not vary much in k 2 at the value m 2 on momentum-squared scales of the order 
of |fc_i_|/a. The contribution of regions of phase space away from this peak can 
be as large or larger as the result due to the peak. In £Q the matrix element 
was split up into a part due to the peak and a part due to the rest of the phase 
space to account for this. Our i|fcjJ/a-prescription gives a good approximation 
of the matrix element away from the peak at k 2 — m 2 , so it makes a more or 
less arbitrary split-up of the cross section unnecesary. 



3.1 Gauge Invariance 

The above prescription breaks gauge invariance. We study the process /i~ + 
/i + — > e~ + D e + W + . To do this, six Feynman graphs with 7, W r± and Z° as 
fundamental bosons are needed. The propagators of the massive bosons must be 
given a width. This does affect the gauge invariance of the amplitude. In [3] it 
was shown that just using the iMF-prescription may lead to grossly inaccurate 
results. However, in this paper we want to focus on the effect of the gauge 
breaking caused by our i|/cjJ/a-prescription. For this reason we use the pole 
scheme for the massive bosons, so that they do not break the gauge. What 
flavour of this scheme we actually used can be found in appendix ^ It turns 
out that in the i?^-gauge, no gauge dependence due to the i|&x|/a-prescription 
is found, although we actually broke gauge invariance. I.e., the results do not 
depend on the gauge parameter £. This can be understood from the Feynman 
graph displayed in equation ^ The gauge dependence comes in via a term 
proportional to (91+92)^(91+92)1/ that occurs in the ^"-propagator. However, 
this term disappears because one of these factors 91 + 92 is to be contracted with 
the current containing the outgoing fermions. These are to be taken massless, 
so consequently this does not contribute, regardless of the gauge breaking that 
may occur at the other side of the VK _ -propagator. To see that our prescription 
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actually breaks gauge invariance we used the axial gauge. In this gauge the 
undressed propagator of the VK-particle is given by 

K z - Mfr + le 

The expression for the squared matrix element can be rewritten in such a way 
that all gauge breaking terms are proportional to | k± | /a or the square of this 
quantity. The axial gauge is not very easy to work with in practice, because 
one either has propagators that mix longitudinal gauge bosons with would- 
be Goldstone bosons or, if propagators are diagonalized, rather complicated 
expressions for the vertices. Details are discussed in _7 4 . To find the gauge 
breaking terms in the unitary gauge we calculate the difference 

\M\ 

gauge-break I I unitary gauge 

-\M\ 

gauge invariant- (11) 

The gauge invariant quantity is calculated by using the axial gauge and the 
gauge invariant prescription 

Res M 

M gauge invariant — , n 9 , ■ 1 7 1 / 

-+M regular V J 

£r — m z + i\k±\/a 

that gives a gauge invariant quantity in the spirit of the pole scheme. This 
calculation was done in the axial gauge to check that we actually get a quan- 
tity that does not depend on the gauge vector n. The algebra was done us- 
ing the C++ computer algebra library GiNaC. C.f., .8. We find that the 
quantity |-M|g auge _ break is, compared to the rest of the cross section, a fac- 
tor |fcjj/(as) ~ l/(ay/s) smaller. Numerically that is a factor 7 • 1CP 14 for 
^/s = 150 GeV and a = ^tt ■ 10 fiva (which is a reasonable value). In ref |3] it 
was shown that gauge breaking effects can get enhanced by a factor as large 
as s/Wg, but even if this would happen, the gauge breaking due to our handling 
of the linear beam size effect remains negligible (note that in the context of 
muon colliders one would actually expect a factor s/m 2 ^ for the case discussed 
in 0). 



4 Conclusions 

The linear beam size effect can be incorporated in the usual Monte Carlo inte- 
gration procedure by doing to substitution 



k 2 — m 2 ± ie k 2 — m 2 ± i\k±\/a 

in the propagator that causes the divergence. This can be done in a gauge 
invariant way, but in the unitary gauge the gauge breaking effect is so small 
that it is safe not to worry about the gauge dependence. The gauge breaking 
effect of the iMr-prescription is much larger than that due to the i\k±\/a. 

Acknowledgement: We would like to thank prof. P. van Nieuwenhuizcn for 
bringing the problem of gauge dependence to our attention. 
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A The Pole Scheme 



To describe resonances, as observed from the W and Z particles, one needs a 
resummed propagator that contains a factor (p 2 — M 2 + iMT) . The prob- 
lem with this propagator is that it breaks gauge invariance, which means that 
observable quantities depend on the gauge choice. The pole scheme (c.f., [S] 
and 0) is one of the ways to solve this. To use it, we first observe that both 
the position of the pole and its residue are gauge invariant quantities. They 
must be because they can be determined by experiment. The consequence is 
that every matrix element that involves such a pole can be written as 

u _ HP 2 = M 2 - iMT) F{p 2 ) - F(p 2 = M 2 ~ iMT) 
J p 2 -M 2 + iMT p 2 -M 2 +iMT +^««ti ^) 

where the first term is gauge invariant, as are the second and third together. 

In practice things are a bit more involved than sketched in the previous 
paragraph. A matrix element generally depends on more that just p 2 and thus 
a prescription is needed to tell us what happens to all the other quantities that 
occur in the matrix element if we put p 2 equal to M 2 — iMT. We follow the 
method outlined in [5] . Our matrix element contains strings of gamma matrices 
with spinors at the beginning and end. These arc canonicalized to ensure that 
all strings of gamma matrices are linearly independent. This means that if we 
have, say, a $ and a ^ in some string of gamma matrices, we can also have the 
same string of gamma matrices with the j> and $ interchanged. We then have 
to decide which of these two comes in front. The anti-commutation relations 
that one has for gamma matrices are then used to do this. Also the relations 
i>u(p) = mu(p) and $v(p) = —mv(p) are used whenever applicable. 

After this has been done, the strings of gamma matrices that remain are 
linearly independent. They are said to form a set of independent covariants. 
If the matrix element is going to be gauge invariant, each coefficient of such a 
string of gamma matrices must separately be gauge invariant. So equation 1141 
is not used for the full matrix element but actually for the invariant coefficient 
that occur in front of the different products of strings of gamma matrices. In 
order to do this, it is also necessary to eliminate one of the outgoing/incoming 
momenta by using momentum conservation for the entire matrix element. All 
inner products between in- or outgoing momenta in the matrix element are 
expressed in a smallest complete set of lorentz invariant variables. In the case 
of the outgoing momenta shown in the graph in equation ^ the set consisting of 

s = (pi +p 2 ) 2 ; 
t= (pi-qi- 92) 2 ; 

x=(qi+q 2 ) 2 ; (15) 

y = (pi +P2 - 92) 2 ; 
z = (pi - qi) 2 , 

can be chosen. If one uses that the squares of incoming and outgoing momenta 
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are given by the appropriate masses squared, all inner products between mo- 
menta are determined by specifying the variables (s, t, x, y, z). Now setting the 
square of some internal momentum in some Feynman graph equal to some value 
is a well-defined operation, except for some caveats that we discuss next. 
The caveats are 

1. If we have outgoing or incoming vector bosons, we should also treat inner 
products of the form p-e with e the polarization vector as linearly independent 
covariant quantities. Some elements in the set of independent covariants 
contain a factor p-e. In the case of the axial gauge this set furthermore 
includes factors that are inner products with the gauge vector n. 

2. In the unitary gauge, the inner product of a polarization vector with the mo- 
mentum of the particle to which the polarization vector belongs is zero. For 
this reason, these inner products should not appear in the set of independent 
covariants, nor in the coefficients that multiply them. The same applies to 
the inner product of polarization vectors with the gauge vector n in the axial 
gauge. 

3. In the axial gauge the property holds that if we have outgoing or incoming 
vector bosons the matrix element becomes zero if the polarization vector of a 
vector boson is replaced by its momentum. It is a feature of the axial gauge 
that this not only happens for massless gauge bosons but also for massive 
ones. This shows that the set of covariants that we had is not really linearly 
independent. To sec how this can be solved consider a matrix element of the 
form 

M = e-piFi + e-p 2 F 2 + e-p 3 F 3 . (16) 

Here the inner products e • pi are the covariants and the Fi are the invariant 
functions. If we know that the relation 

q-p 1 F 1 +q-p 2 F 2 + q-p 3 F 3 = (17) 

holds, we can eliminate F\ from the matrix element. We get 

M= (e-p 2 -e- Pl ^-)F 2 + ( e .p 3 -e-Pi— W ( 18 ) 

V q-pij V q-pij 

Thus we have actually reduced the set of covariants from three to two in this 
example. This boils down to doing the substitution 

E^e — E-pi — - — . (19) 

q-pi 

In this substitution the vector p\ can be chosen to be any linear combination 
of incoming and outgoing momenta. It is advisable to choose one that does 
not yield any singularities in the physical phase space due to dividing by q-pi- 
In the unitary gauge a similar reduction can be carried out. In our calculation 
we chose to get (qi + q 2 ) ■ q 3 in the denominator. This is equal to s — x — . 
This quantity has no poles in the physical phase space. 
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4. One has to be careful about the set of invariant variables. Actually the 
set (s, t, x, y, z) has a problem. To see this, consider the Feynman graph 

(20) 




M (PlY e (gi) 

The internal electron propagator is given by 

s = _Mi±M_. (21) 

s-x-y + M^ 

In the pole scheme we should determine the residue for the Z-pole. This 
means that we put s = M| to lowest order. The maximum value of x + y 
is s and occurs in the limit that the outgoing electron is produced at rest. 
We see that the quantity l/(s — x — y + M^) does not have a pole in the 
physical phase space but if we put s = M§ it does develop a pole. For 
this reason we did not use the set of parameters (s, t, x, y, z) but instead the 
set (s,t, £, rj,z) where £ = x/s and 77 = y/s. This set causes no trouble with 
spurious singularities. 
The problem with spurious singularities, that we encountered in item 3 and 4 
can be looked upon as follows. Formula^]tells us to split the invariant functions 
in the matrix element. However, we have some freedom in making this split- 
up. This makes it possible that the pole term has a singularity that is then 
canceled by the regular term. A sensible choice of such a split-up takes care not 
to introduce new singularities in the physical phase space. 
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